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ABSTRACT 

We present IRAM Plateau de Bure interferometric detections of CO (/=! — s-0) emission from a 24/im-selected 
sample of star-forming galaxies at z = 0.4. The galaxies have PAH7.7/im-derived star formation rates of 
SFR^30-60Mo yr"' and stellar masses M* ^ lO^M©. The CO (7=1— ^0) luminosities of the galaxies im- 
ply that the disks still contain a large reservoir of molecular gas, contributing ^20% of the baryonic mass, 
but have star-formation 'efficiencies' similar to local quiescent disks and gas-dominated disks at z L5- 
2. We reveal evidence that the average molecular gas fraction has undergone strong evolution since z ~ 2, 
with /gas oc (1 The evolution of /gas encodes fundamental information about the relative deple- 

tion/replenishment of molecular fuel in galaxies, and is expected to be a strong function of halo mass. We show 
that the latest predictions for the evolution of the molecular gas fraction in semi-analytic models of galaxy 
formation within a ACDM Universe are supported by these new observations. 

Subject headings: galaxies: evolution — cosmology: observations 



1. INTRODUCTION 

Molecular hydrogen is arguably the most important compo- 
nent of the interstellar medium (ISM) in high-redshift galax- 
ies, since it is the phase necessary for, and immediately pre- 
ceding, star formation. However, observing cool (<50K) 
molecular hydrogen is difficult; the molecule has no perma- 
nent electric dipole and so one must rely on tracer molecu- 
lar emission (the most common being ro-vibrational emission 
from the '^CO isotopomer [hereafter 'CO']). 

Unfortunately, uncertainty of the precise calibration of CO 
luminosity to total gas mass in both low- and high-z galax- 
ies have made interpretation of results challenging; even in 
the local Universe there is seen to be metallicity and lumi- 
nosity dependence on the conversion (Boissier et al. 2003; 
Blitz et al. 2007; Komugi et al. 2010). Furthermore, although 
the bulk of the molecular gas in star-forming disks is con- 
tained in virialised clouds, it is clear that the most luminous 
galaxies - mostly mergers - have molecular gas distributions 
and star-formation modes dramatically different to quiescent 
disks; they require an alternative molecular mass calibration 
related to the physical conditions of their interstellar media 
(Solomon et al. 1997). The importance of understanding how 
the thermodynamic state of the gas affects the calibration 
also remains controversial (high-z observations typically tar- 
get /upper > L requiring some correction to 7 = (1 — > 0) based 
on the excitation of the gas; although see Ivison et al. 2010; 
Harris et al. 2010). 

The majority of molecular gas studies at high-z have, nec- 
essarily, focused on the most active systems (generally sub- 
millimeter selected galaxies and quasars), and so compara- 
tively little is known about the evolution of the molecular gas 
properties of the more common, but less active population. 
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Nevertheless, recent studies have started to make progress in 
expanding the parameter space of CO observations at high-z, 
and several studies have detected CO emission from 'normal' 
disks at z = 1-2 (Daddi et al. 2008; 2010; Tacconi et al. 2010). 
The striking result it that the efficiency of star-formation in 
these distant galaxies is claimed to be similar to that seen 
in local 'quiescent' disks (Dannerbauer et al. 2009; Daddi et 
al. 2010), with enhanced SFRs driven simply by larger gas 
reservoirs, drained by star-formation (and feedback) on Gyr 
timescales. This hints that a secular mode of star formation 
might be ubiquitous in L* disks independent of cosmic epoch 
(but see the higher-z results of Coppin et al. 2007 and Riechers 
et al. 2010). 

The aim of this Letter is to further expand the parameter 
space of molecular gas studies in normal galaxies, building 
on our observations of luminous infrared galaxies (LIRGs) 
at z = 0.4 (Geach et al. 2009, G09). We present new IRAM 
Plateau de Bure interferometric (PdBI) observations of a sam- 
ple of 24/xm-selected galaxies in the outskirts of the rich clus- 
ter CI 0024-1-16 at z = 0.395. These galaxies bridge the gap 
between quiescent, local spirals and luminous, high-redshift 
starbursts and disks. Throughout we adopt a 17^ = 0.3, f^A = 
0.7, Ho = 70km s"' Mpc"' cosmology. 

2. OBSERVATIONS AND DATA REDUCTION 

Our sample consists of five Spitzer MIPS 24/im detected 
LIRG-class galaxies in the outskirts of the rich cluster 
C10024H-16 (z = 0.395, see Geach et al. 2006). We obtained 
mid-infrared spectroscopy with Spitzer IRS to confirm the 
presence of aromatic features, a clear indication of star for- 
mation (rather than infrared emission dominated by an active 
galactic nucleus; Geach et al. 2009b). Combined with the 
two galaxies observed with IRAM PdBI in G09, this repre- 
sents an effort to obtain CO detections or limits for all the 
galaxies where we have detected significant emission from 
the polycyclic aromatic hydrocarbon 7.7/im band. The line 
luminosities imply star formation rates of ~3O-6OM0 yr"', 
thus allowing us to link molecular gas with moderate levels of 
dusty star formation in the galaxies. 

In this Letter we present the results in the context of the 
broader star-forming population at z = 0.4, rather than focus- 
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Table 1 

Properties of the galaxy sample 



Target 


"J2000 


<5j2000 


Z[0 II] 


^FIR 


SFR 


"^FWHM 


^CO(l-O) 


M(H2+He)t 


mV 


/gas 




(h m s) 


(o ' ") 




(1O"L0) 


(Meyr-') 


(kms-') 


(10'" K km pc^) 


(IO'OMq) 


(10"Mq) 




MIPS 1002652.5 


00 26 52.5 


+17 13 59.9 


0.3799 


3.5±0.5 


62 ±19 


140 ±10 


0.64 ± 0.05 


2.9 ±0.2 


0.95 


0.23 ± 0.02 


MIPS J002703.6 


00 27 03.6 


+17 11 27.9 


0.3956 


2.3 ±0.3 


42 ±12 


250 ±30 


0.43 ±0.06 


2.0 ±0.3 


0.87 


0.19 ±0.03 


MIPS J002715.0 


00 27 15.0 


+17 12 45.6 


0.3813 


1.9 ±0.3 


35 ±11 


340 ±40 


0.26 ±0.03 


1.2±0.3 


1.10 


0.10 ±0.03 


MIPS J002609.1 


00 26 09.1 


+17 15 11.5 


0.3940 


1.5 ±0.3 


28±9 




<0.39 


< 1.8 


0.41 


<0.30 


MIPS J002606.1 


00 26 06.1 


+17 04 16.4 


0.3904 


1.9 ±0.3 


34 ±10 




<0.25 


< 1.1 


0.49 


<0.19 


MIPS J002621.7* 


00 26 21.7 


+17 19 26.4 


0.3803 


3.1±0.2 


56± 16 


140 ± 10 


0.68 ±0.06 


3.1 ±0.3 


1.12 


0.22 ± 0.02 


MIPS J002721.0* 


00 27 21.1 


+16 59 49.9 


0.3964 


3.2 ±0.2 


59± 16 


160 ±30 


1.14±0.11 


5.2 ±0.5 


0.98 


0.34 ±0.04 



I Assuming a CO luminosity to total gas mass conversion factor of a = 4.6M0(Kkm pc ) that includes helium 

tt Uncertainties '^0.02 dex based on refitting spectral templates 1000 times with random perturbation within \a photometric uncertainties. 

*Geachetal. (2009) 



ing on potential environmental effects in the cluster. We be- 
lieve this is a valid approach, since these galaxies are expected 
to be randomly accreted onto the cluster from the surrounding 
field, and the galaxies are seen at a stage where strong envi- 
ronmental effects associated with clusters (specifically ram- 
pressure stripping and harassment) are yet to have an effect 
(Treu et al. 2003; Moran et al. 2007). They lie at clustocentric 
radii of 1 .8^.6Mpc (~ 1-3 x the virial radius of CI 0024+16). 
Intrinsically, they are likely to be representative of disk galax- 
ies of equivalent mass sampled from a random field, and can 
be compared to galaxies with similar properties at low and 
high redshift; i.e. normal star forming disks. 

The new IRAM observations were conducted over June- 
October 2009 in configuration 'D', using 5 antennae, and we 
adopted the same strategy as in G09. Sensitivities ranged be- 
tween 0.5 1-1 . 1 7 mJy (median average for 10 MHz wide chan- 
nels and 2 polarisations), and the on-source exposure times 
ranged between 6.4-12.8hrs. We targeted the CO(7=1->0) 
115.27GHz rotational transition at ;/obs — 82.63 GHz in 5 
sources. The 3 nam receiver was tuned to the frequency of 
the redshifted CO (7=1^0) Une at the systemic redshift of 
each galaxy derived from optical spectroscopy (Czoske et 
al. 2000; Moran et al. 2007). As in G09, the correlator 
was set-up with 2.5 MHz spacing (2x64 channels, 320 MHz 
bandwidth). The phase and flux calibrators were the sources 
3C454.3, 0119+1 15 and 0007+171. The observing conditions 
were good or exceUent in terms of atmospheric phase stability, 
however any anomalous and high phase-noise visibiUties were 
flagged in the calibration stage. Data were calibrated, mapped 
and analysed using gildas (Guilloteau & Lucas 2000). 

3. RESULTS 

3.1. Plateau de Bure CO detections 

We detect CO (7=1^0) emission in three galaxies (>4(t de- 
tections within 2" of the phase tracking centre, significances 
determined from the integrated line flux). Velocity-integrated 
maps and nnm-spectra are shown in Figure 1. Total fluxes are 
evaluated from single or double Gaussian fits to the spectta 
depending on whether the profiles are very broad and have 
hints that they are double peaked. The line luminosities are in 
the range L'^q = 0.26-0.64 x 10"' K km s'^ pc^ (Table 1). Un- 
certainties on the luminosities are estimated by re-evalutating 
the Gaussian fits after artificiaUy adding noise to each chan- 
nel based on the observed r.m.s. fluctuations in each data cube. 
Upper 3(T limits on L(^q for the two non-detections are based 
on the r.m.s. noise in the observations of each source (note 
that the two non-detections are also the two galaxies with the 
lowest SFRs in the sample). The luminosities and Une widths 



are listed in Table 1; we also list the properties of the two 
galaxies from G09 which are included in the following analy- 

sis 

The spectra of MIPS J002652.5 and MIPS J0027 15.0, Uke 
MIPS J002721.0 (G09) exhibit broad CO (7=1^0) emission, 
and can be well-fit by double gaussian profiles. None of the 
new sample show obvious signs of major mergers or strong 
tidal interaction in the deep optical imaging (Fig. 1), and so 
we conclude that we are most likely observing CO (7=1^-0) 
emission tracing molecular gas distributed over rotationally- 
supported disks. 

3.2. Physical properties of the galaxies 

Stellar masses for aU the galaxies were determined by 
fitting the galaxy SED constructed from BVRIJK imaging 
(Moran et al. 2007) to a large suite of model SEDs using the 
KCORRECT software, v4.2 (Blanton et al. 2007). Utilizing 
the known spectroscopic redshifts, KcORRECT finds the best 
non-negative combination of the template spectra to fit the 
galaxy SED in its rest-frame, including a reddeiung law, with 
best-fit values covering a range Ay = 0.8 -2.3 mag. Stellar 
masses are then calculated from the luminosity and mass-to- 
light ratio of the best-fitting model in the rest-frame K-band, 
= Mq/Lq. The fits had an average mass-to-light ratio 
of (T*) = 0.36 (range 0.24-0.52), and inferred stellar masses 
ranging overM^ ~0.5-l x IO'^Mq (Table 1). 

In G09 we adopted a conservative 'ULIRG' conver- 
sion of a = 0.8MQ(Kkms"' pc^)"' for gas mass, where 
M(H2 + He) = uLq (we omit the units of a for clarity in the 
following). However, there are hints that a Galactic scaUng is 
more appropriate for these galaxies. The first point to note 
is the similarity between the so-called star-formation 'effi- 
ciencies', measured by L'^q/Li^, compared to local quiescent 
disks. The disks in this sample have a mean (Lir/L(,q) =51 
(range 28-73), similar to that of local spirals, but lower than 
that of ULIRGs, which tend to have T-ir/T-co ^100 (Solomon 
et al. 1997). Daddi et al. (2010) and Tacconi et al. (2010) re- 
port similar 'secular' properties for z > 1 disks, although note 
that these high-z observations measure 7up > 1 transitions, and 
so require an additional uncertain correction to estimate gas 
mass. It is claimed that the efficiency of star-formation re- 
mains relatively constant in typical disks at all epochs (al- 
though, as a concept, star-formation 'efficiency' only makes 
sense if all of the molecular gas has an equal probability of 
undergoing star-formation). 

Although the CO (and mid-infrared) emission is unre- 
solved, we can evaluate how appropriate a = 0.8 is based 
on the offset from the Keimicutt-Schnudt (K-S) law (Ssfr oc 
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Figure 1. Velocity integrated, clean CO (7 =1—^0) line emission maps (at levels of >2 X r.m.s. in steps of 0.5) overlaid on 30" X 30" Subaru SuprimeCam 
i?-band images (including a 5" X 5" HST/WFPC2 inset for MIPS J002703.6). The ellipse indicates the beam size and shape, (top) 3 mm spectra extracted from 
the peak of the integrated line emission maps, fit with either a single or double Gaussian (where doubles have a fixed amplitude and width for each peak). 

a mean of (/gas) = 0.21 ±0.03. We can conclude that these 
are actively star-forming galaxies with large gas reservoirs, 
but the disks are dominated by stellar mass. This is markedly 
different to the recently discovered z ^ 1-2 gas-rich disks of 
Daddi et al. (2010) and Tacconi et al. (2010), which are gener- 
ally more active galaxies, with comparable stellar masses and 



Sgal; Kennicutt 1998), assuming both gas and star-formation 
trace the stellar emission. Inclination and effective radii of the 
disks are estimated by fitting a 2D model following a Sersic 
profile to a CFHT /?-band image (Czoske et al. 2001), using 
GALFIT v3.0.2 (Peng et al. 2010). GALFIT was run on postage 
stamps of 20" x 20", simultaneously fitting models to the 
main galaxy and any bright companions in the field. GALFIT 
produces estimates of the galaxy axial ratio and size that ac- 
counts for the seeing of the ground based image, by convolv- 
ing with a point-spread function that we have extracted from 
stars in the CFHT image. The effective radii of the disks range 
between 0.9-1.5" (4.8-8. Ikpc), with a mean of 1" (5.4kpc), 
and typical uncertainties of 0.02". We take the effective radii 
as the size of the disks and find a mean (q:ks) ~ 4.2, in line 
with a Galactic conversion. Note that this should be taken 
only as a guide for choosing between a Galactic and ULIRG- 
like conversion, oks should not be applied as the actual cali- 
bration: it is particularly sensitive to the assumed spatial dis- 
tribution of the gas and star-formation, such that a larger as- 
sumed radius would result in a larger oks and vice versa (with 
r < 1 kpc required for a^s ~ 0.8). 

MIPS J002703.6, is partially covered by one of the sparse 
Hubble Space Telescope/WFPCl (F814W) images taken of 
CI 0024H-16 (Treu et al. 2003), reveahng finer detail of the op- 
tical morphology (Figure 1). The galaxy has been typed as 
an early-type spiral (Treu et al. 2003), and has clear spiral 
structure with several bright knots of emission in the northern 
arm (presumably star-forming regions) and a potential bar. 
We have run GALFIT on the HST/WFPC2 image, where the 
best 2D fit is a two-component Sersic bulge and exponential 
disk, where the bulge and disk have effective radii of 0.15" 
(0.8 kpc) and 1 .05" (5.6 kpc) respectively (consistent with the 
value derived from the ground-based imaging), and the bulge- 
to-total Hght ratio is B/r = 0.03 1 . 

Based on this evidence, we now assume a Galactic con- 
version of a = 4.6 to estimate the gas masses, which implies 
molecular gas masses of M(H2 + He) =1.2-2.9 x lO'^M©. 
What is the gas contribution to the total baryonic mass of 
these galaxies? With our choice of a, the resulting gas frac- 
tions, /gas = Mgas/(Mgas+M*), are in the range 0. 10-0.34, with 



higher gas-fractions /gas > 40% 

4. INTERPRETATION 

Figure 2 shows the observed /gas in > 1O'°M0 galax- 
ies at z = (Leroy et al. 2008), z = 0.4 (this work) and 
z= 1.5-2 (Daddi et al. 2010, Tacconi et al. 2010) showing 
the clear decline in (/gas) in the lOGyr since z ^ 2. As a 
guide, we cast the evolution of /gas relative to the cosmo- 
logical abundance of baryons in stars and molecular gas at 

z = 0: f^H, /(^K2 + ^*) = O-'^^^-oim from the 'baryon census' 
of Fukugita, Hogan & Peebles (1998). At z = 0, the observed 
/gas in star-forming disks is close to that inferred from the total 
relative abundance of molecular gas and stars in the Universe, 
but at z 2 it was a factor ^5-10 larger than the local valueQ. 
The decline in molecular gas fraction in the stellar mass lim- 
ited sample is broadly characterised by /gas oc (1 +z)'', with 
7^1 .5-2.5, shallower than the rate of decline of the average 
specific SFR in galaxies over the same epoch, which falls off 
more like (1+z)^ (e.g. Karim et al. 2010). This is expected, 
since molecular reservoirs can be replenished with fresh mate- 
rial via cooling; thus, it is the relative evolution of the cooling 
rate and star formation rate that shape the evolution of /gas. 

Recent semi-analytic prescriptions for galaxy formation 
make predictions for the evolution of the molecular ISM. In 
Fig. 2 we show the average Hi + He fraction in > 10^^ Mq 
galaxies selected in halos of various mass from the GALFORM 
model (galaxies populated within the Millennium Simulation 
ACDM framework). This latest model (Lagos et al. 2010) 
implements an empirically based star-formation law (Blitz & 

^ The similarity of the average galaxy gas fraction and the cosmic abun- 
dance at z = is perhaps coincidental, since is integrated over the lumi- 
nosity functions of all bulges, disks and iiTegulars, whereas is inferred 
from a relatively limited local CO survey. However, it could be a convenient 
datum for comparisons of /gas. 
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Figure 2. Evolution of the molecular gas fraction in > IO'^Mq galaxies, showing decline in /gas over cosmic time following (1 +2)^^^"^. We compare our 
z = 0.4 results with the high-z disks of Daddi et al. (2010) and Tacconi et al. (2010), and the z ~ sample of Leroy et al. (2008). Large open symbols show the 
average of each sample, with error-bars indicating the range. The upper bar of the Daddi et al. sample shows the effect of assuming a = 4.6. All stellar mass 
estimates assume a Chabrier initial mass function (IMF; we have applied a factor 1.25 to the Leroy et al. data). We also show the evolution of /gas in galaxies 
with M* > IO'^Mq in different halos from the latest semi-analytic models (Lagos et al. 2010). 



Rosolowsky 2006) to estimate the molecular gas mass. The 
model predicts that at z < 3 /gas is lower and the rate of de- 
cline faster for more massive halos, and at in the mass regime 
pertinent to the samples presented, 1O'^M0, evolves in a way 
similar to the observations. The results are broadly consistent 
with other models. For example, in smoothed particle hy- 
drodynamic simulations, Keres et al. (2005) predict that the 
average accretion rate of gas onto galaxies within halos of 
Mvii- = IO'^Mq decreases from Mcoid ^ 20 Mq yr"' at z = 2 to 
Mcoid ^ 3 Mq yr"' at z = 0.4. The observed evolution is also 
consistent with the predicted cosmological evolution of the 
accretion rate at fixed halo mass in other models (see Dutton 
et al. 2010). 

The evolution of /gas is expected to be strongly dependent 
on galaxy mass, since the evolution of the global average SFR 
appears to be halo mass dependent, as does the expected evo- 
lution of the gas accretion rate. Cooling and star formation 
are the key drivers of the molecular gas fraction; but it is 
also sensitive to the merger rate (which can deliver molecular 
material, trigger star-formation and reconfigure the baryonic 
content), feedback (heating or ejection of cold gas and H2 dis- 
sociation), environmental effects (e.g. ram-pressure stripping) 
and growth in size of galaxy disks. The latter has an impor- 
tant effect since it (on average) reduces the disk hydrostatic 
pressure (oc r^) and therefore the efficiency of H2 formation 
(reflected in the H2/H1 ratio, Elmegreen 1993; Wong & Blitz 
2002), as can the metallicity and strength of the interstellar 
UV radiation field. 

Selection effects will affect Fig. 2 severely, and we should 



discuss how these might affect our interpretation. Sensitivity 
limits do not allow us to detect galaxies with very low gas 
fractions at high-z, and in low-z surveys the very gas rich sys- 
tems appear to be rare. Here we have attempted to compare 
similar galaxies across a range of epochs, but the current sam- 
ple sizes are small, and inhomogeneously selected; any con- 
clusions we draw from Fig. 2 should be treated as indicative 
for future surveys that can more precisely control selection 
and have the efficiency to survey larger numbers. 

A further caveat to note is that both Daddi et al. (2010) 
and Tacconi et al. (2010) apply corrections to their ob- 
served CO luminosities on account of the difference in the 
Rayleigh-Jeans brightness temperature of the higher J transi- 
tions targeted. Tacconi et al. (2010) apply Rt,i = 0.5 to their 
CO (7=3— !>2) luminosities (defined Rt,[ = ico(3^2)/^co(i^O))' 
and Daddi et al. (2010) apply R21 = 0.6 to their CO (7=2-^1) 
luminosities (see also Dannerbauer et al. 2009). Leroy et al. 
(2008) also apply R21 = 0.8 to a sub-sample of the z ~ galax- 
ies that were mapped in CO (7=2^-1). Furthermore, Daddi et 
al. (2010) assume an a = 3.6 ±0.8, inferred from dynamical 
mass arguments. Although this is within ~1ct of the Galactic 
value, assuming a = 4.6 increases the gas fraction by ^12%. 
These systematic uncertainties in CO-H2 conversions should 
be taken as a corollary when interpreting Fig. 2, however a 
large change in the gas mass (or stellar mass) would be re- 
quired to remove the trend in /gas. It is unlikely that we have 
over-estimated the gas mass in the z = 0.4 LIRGs, since we 
have applied a = 4.6. In the high-z sources, the gas masses 
would have to be over-estimated by a factor ^^5, or the stellar 
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masses under-estimated by a factor ~4 to remove the trend. 

5. SUMMARY 

We have presented new constraints on the molecular gas 
mass of z = 0.4 LIRGs, showing strong evolution of /g^s since 
Z^2. The observed trend encodes critical information on the 
evolution of the growth of galaxies, including gas accretion, 
feedback and stellar mass assembly. Future surveys that are 
able to split samples based on environment, mass, and other 
parameters will disentangle the relative effects that shape the 
evolution of the gas fraction. These observations will soon be 
supplemented and surpassed by progress made with the Ata- 
cama Large Millimeter Array and the Expanded Very Large 
Array, which will provide unprecedented sensitivity and reso- 
lution with which to probe the evolution of gas in galaxies. 
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Sciences and Engineering Research Council of Canada, sup- 
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trophysics and Cosmology (McGill) and the U.K. Science and 
Technology Facilities Council (STFC). I.R.S. and A.C.E also 
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